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1. Introduction 

The investigation of exact solutions to nonlinear evolution equations plays an important 
role in the study of nonlinear physical phenomena. The wave phenomena are observed 
in physics, mechanics, biology, etc. The effort to find these solutions is significant for 
the understanding of many physical phenomena, thus they may give more insight into 
the physical aspects of the problems. 

This paper is devoted to study the solutions of the variable coefficient Gardner 
equation which is given by 

U t + HiUU x + [JL2U 2 U X + flsU xxx + h(t)u x = 0, (1) 

where /ii, /i2 and ^ are constants and h{t) is a function of t. Equation (CQ) includes 
considerably interesting equations, such as KdV equation, when /i 2 and h{t) are equal 
zero, and mKdV equation, when {i\ and h{t) are equal zero. In arterial mechanics, some 
special cases of Equation ([!]) were considered as follows: 



(i) The variable coefficient KdV equation 

u t + ii\uu x + (J, 3 u xxx + h(t)u x = 0, (2) 

was considered in pHF]. 

(ii) The variable coefficient modified KdV equation 

u t + ^u 2 u x + fisu xxx + h(t)u x = 0, (3) 

was considered in [HIT]. 

(iii) The KdV equation 

u t + [i\uu x + [i 3 u xxx = 0, (4) 

was considered in [8HT9] . 

(iv) The modified KdV equation 

u t + ^u 2 u x + [i 3 U xxx = 0, (5) 
was considered in 118-201. 



Special cases of have been studied in [21TT25] . The rest of this paper is arranged 
as follows. In section [2j we apply the Lie classical symmetry analysis method to (CE]). In 
Sec. [31 the mathematical framework of the generalized Jacobi elliptic function expansion 
method will be provided. In Sec. HI some new exact solutions of (CQ) are obtained and 
application of the obtained solutions in arterial mechanics will be presented in section 

El 
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2. Classical symmetries 



To apply the Lie classical symmetry analysis method [261ET] to ([I]), we consider the 
one-parameter Lie group of infinitesimal transformations in x, t, and u given by 

x* =x + e£(x,t,u) + 0(e 2 ), (6) 

t* = t + er(x,t,u) + 0(e 2 ), (7) 

u* = u + er)(x,t,u) + 0(e 2 ), (8) 

where e is the group parameter. We require that the set of solutions of (CQ) be invariant 

under this transformation. This yields an overdetermined system of linear equations 

for the infinitesimals £(x,t,u), r(x,t,u), and r)(x,t,u). The associated Lie algebra of 

infinitesimal symmetries is the set of infinitesimal generators of the form 

Odd 
v = £(x, t, u)— + t(x, t,u)— + r](x, t, u))-^- (9) 

Invariance of ([1]) under a Lie group of point transformations with infinitesimal generator 
()9]) leads to a set of determining equations. Solving this system we will obtain: 

r = Qciii 2 t + c 2 , (10) 
f = cl \2^ 2 x - n\t + 6[x 2 th(t) - 2/x 2 J h{t)dtj + c2h(t) + c3, (11) 

77 = -cl(2^ 2 « + /i 1 ), (12) 

where, ci, c 2 and C3 are arbitrary constants. The associated infinitesimal generators are 
given by: 

(r \ 

2fx 2 x - jj^t + §n 2 th(t) - 2ji 2 J h(t)dt\ — + 6^2^ - {2/j, 2 u + (13) 

v 2 = A, (14) 

«*-*We + I- (15) 

When considering the infinitesimal generator i>!, we will obtain the surface condition 

2fi 2 x - n\t + 6fi 2 th(t) - 2/j 2 J h(t)dt) fa+S^t-fo = -(2/i 2 ^+/Ui), (16) 

which when solving we will obtain the similarity transformation 

2 

u = ^ + t^f(C), ( = xt^ + p-tl - / /l(t)tft, (17) 

when substituting from (TTT1) into ([1]) we will obtain the following ordinary differential 
equation 

-/-C/' + 3/X2/V' + 3//3/'" = 0. (18) 

In general, the exact solution for this nonlinear ODE can't be obtained by using the 
elementary functions. The series solution of this equation was obtained in [28J. 
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When considering the infinitesimal generator v 2 + ctf 3, a is a constant, we will obtain 
the surface condition 

Wt)+ a)| + | = 0, (19) 
which when solving we will obtain the similarity transformation 

u = f(QX = x- J h(t)dt-at, (20) 

when substituting from (1201) into ([1]) we will obtain the following ordinary differential 
equation 

- af + frff + fi 2 f 2 f + Ms/'" = 0, (21) 
this equation will be solved in the following section. 

3. Generalized Jacobi elliptic function expansion method 

For solving (1211) we will use the generalized Jacobi elliptic function expansion method 
[29l[30]. It is assumed that (121]) has the solutions in the form 



u 



a + E a i<P\ ( 22 ) 



where a, are constants to be determined later and <f> satisfy the following elliptic 
equations 

cf)' 2 = r + p ^ + q( p\ (23) 
furthermore we can get 

0"=P0 + # 3 , (24) 

where the primes denotes the derivatives with respect to £ and r, p, q are constants. 
The solutions of (1231) are written in Appendix A [31J. Balancing the highest derivative 
term /"' with nonlinear term //' in (j2ip gives n — 1, from which we have 

u = ao + a\(f) + a_i0 _1 . (25) 

Substituting (|25|) along with (123!) and (124"]) into (I2T]) and collect all terms with the same 
powers in <j) k (p n (k = 0, 1, 2, ...,/ = 0, 1, ...) and set the coefficients of these terms to zero 
yields a system of nonlinear algebraic equations for a , a± and a_i. Solving this system, 
we find three sets of solutions. 
The first set: 



a = - — , ai = 0, a_i = ±W , a = . (26) 

2/i 2 V ^2 4// 2 

The second set: 



a = - — ,01 = ±4/ ,a_ 1 = 0,a = . (27) 

The third set: 



-/ii . l-6qU3 /-6ru 3 -4/x 2 /i3(±6v/rg - p) - //? 
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4. Exact solutions 



In this section we will give some exact solutions for (JTJ. Substituting from ( 1261) . ( 1271) and 
( 128]) into fl25|) and using the table in Appendix A we will obtain the following solutions 
of©: 

Case 1: Soliton and soliton-like solutions 



-^i^H^M^W- (29) 



B2 = Zffl ± & ech ^ _ / h{t)dt _ f ^) A . (30) 

2/i 2 V Ma V ^ V 4 A*2 / / 



W3 = ^ ± ,F^tanh (x - [ h(t)dt - ( Z^M^A ) t ) . (31) 



2/i2 V A^2 v y V 4^ 



--S^-H^/^-pS^ > • (32) 



u 5 = — — ± J — — tanh(x- / h(t)dt - at) ± J — — coth(;c - / h(t)dt-at), (33) 

2^2 V ^2 •> V ^2 •/ 

with 

tt= /-w±6+2)-,n (34) 

Case 2: Triangular periodic solutions 



w 6 = 


2^2 V 


/-6/Xa 

fJ-2 


«7 = 


2/i 2 V 


/-6// 3 


u s = 






u 9 = 


2/i 2 V 




u w = 


2^2 


/-6/X3 
V ^2 



sec ( x - / fc(f)<ft - f " 4 ^"^ ) tj . (35) 

esc ( x - / - t \ . ( 36 ) 

tan ( x - I h(t)dt - ( 8 ^~ M * ) f) • (37) 

cot ( x - / fc(i)<ft - f 8 ^"^ ) *) • (38) 



tan(x - J h(t)dt - at) ± ^ -^-cot(x - J h(t)dt - at), (39) 



with 

a = (- 4w( ^; 2) -" f ) . (40) 
Case 3: Jacobi elliptic function solutions and combined Jacobi elliptic function solutions 



— ir * iff - * " / " ( ~ w 4r~" ? ) ' ' (41) 
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„ 13 = ± EM^ nc L _ / ft(t)di _ ( Watfc' ~ n - A (43) 



2^2 V ^2 V ^ V 4 // 



'-2 



-1 + fcsn [ x - / /i(t)cft 



"» = ^ * v — — dn (_ /Mt)(it - - • (47) 



\ 4/^2 



-1 + sn (x-fh(t)dt-( 4 ^7f-* ) t 



u i8 = \rr ± y— a t 77. - 1+fc2 2 x x • ( 48 ) 



-1 + dn (x-fh(t)dt- ( z^t^7A ] t 



/Z^ k' + dn 2 (x - / h(t)dt - ( -^M3(i+fe' 2 )-M? \ A 

2^ 2 V ^2 k> - dn 2 (a - J fe(t)dt - ^ t 



± dn 2 (x - / - ( teO^bg t ) 



M21 "^2" V^T 2VFdn (x - J fe(t)<ft - ( W^EE g)^) • 1 J 



, , 2 cn 2 (x - Jh(t)dt - at) + k'sn 2 (x - J h(t)dt - at) 

M22 ~2^ 2 V ^2 } cn 2 (x - J h(t)dt - at) - k'sn 2 (x - J h(t)dt - at)' 1 j 

with 

a = (-^(^P 2 ~" 1 ) (53) 



/-6^ 3 (l + *0 : 

u 23 — T, ± 



2^2 V ^2 



2 (l ± fan (x - J h(t)dt - ( Zggjj|±g^ t )) 
*\ (1 + fc) (l ± sn (x - J fe(t)dt - ( z^+gl^I ) t )) • 
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J(l + A;sn(:r - / h{t)dt - at))(l + sn(x - / h{t)dt - at)) 

cn(x - J h(t)dt - at) ' v ; 

with 

a = + . (56) 



2^2 V A*2 



^(../^.( W^-rf ),). (59) 



- - ^ ± ^ (• - / ^ - ( 4 " 2 " 3(i :r ) "" f ) <) ■ < 6 °» 



—^ ± ^ to ('-/^-( w y ) ~ ij )')- < 6i » 



1 + fen f x - J h(t)dt - 



«3i = ^±J— ^ S S — ^ — \~ - (63) 



,(l±^!l_ u 2 

cn I x — : ' ■ ' ' ■ • ■ / ' - ' — — — 1 



4^2 

_ '21 



- = ^ ± V— T ; + - r _ - — _ ( - w^-, — ■ (65) 



35 2„ 2 * V ^2 ( + ' k> ± dn 2 (*. - / h(t)dt - ( ^Cf^)--j ) t ) • (t> ° 
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-fJ-i , / -6jU 3 cn 2 (x - Jh(t)dt - at) - k'sn 2 (x - / h(t)dt - at) 

Um ~ 2/i 2 V ^2 cn 2 (x- fh(t)dt- at) + k'sn 2 (x-fh(t)dt -at)' 1 j 

with 

a = (- 8 ^' 1 4 + fc, » 2 -"') . (69) 



-//i /-6^3 



(1 + jfc) (l ± sn (x - J h(t)dt - ( -/WW* 2 )-/*? ) t 
^2 '"J 2 (l ± fen (x - / h(t)dt - ( -/WihW)-/£ ) 



-/il -6^3 / l - fc 



with 

■At2A*3(l -6& + A; 2 ) - /if 



ct = 

M39 = 
with 
a = 

U40 = 
with 

a = 

U41 = 
with 

a = 

U42 = 
with 

a = 



4p 



2 



7/1 ± W-^ns(x - I - at) ± J-^-k 2 sn(x - J h(t)dt - at), (73) 



2/t 2 V A*2 J V A*2 

-4^3(±6v^+l + A; 2 )-Atf 



4/i 



2 



-111 \ I'! J " V /'2 

-4^ 3 (±6v^ + 1 + A; 2 ) - fij 



4/i 2 



(70) 



x cn(x - J h(t)dt - at) 

^(1 + ksn(x - / h(t)dt - at))(l + sn(x — / h(t)dt — at)) 



(72) 



(74) 



(76) 



^ ± JJ^lk' 2 nc(x - J h(t)dt - at) ± J^-k 2 cn(x - J h{t)dt - at), (77) 



2^2 V ^2 ^ V ^2 



-4/i 2 /x 3 (±6v /z F^ + A; ' 2 - k 2 ) — /if \ g 



4/. 



2 



± JJW*cs(x - f h(t)dt - at) ± J—^k' 2 sc{x - f h(t)dt - at), (79) 

2/i2 V ^2 V ^2 J 

-4/i 2 /i 3 (±6v^-l-^ 2 )-/x 2 \ 

4// 2 
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u 43 = — ± J^k' 2 nd(x - f h(t)dt - at) ± J^dn(x - f h(t)dt - at), (81) 

2/i2 V ^2 J V ^2 J 

with 

a = / -^Mie^-i-^-^ j ^ (g2) 



M 44 = — ± J—k 2 k' 2 sd(x - I h(t)dt - at) ± J-^-ds(x - / h(t)dt - at), (83) 

2/i2 V ^2 V ^2 J 

with 

/ 4^(±6v^Fg - 1 + 2A:' 2 ) - /i 2 ^ 



—fii — 6/j 3 k 2 — 1 1 + &sn(:r — / h(t)dt — at) 
M45 2/i 2 V ^2 4 dn(x - / h(t)dt - at) 

I — 6/13 A; 2 — 1 dn(a; — / h(t)dt — at) 
\]~ 2 4 1 + A;sn(x - / /i(t)dt - at) ' ( ^ 

with 

/-w±f- |)y-D-^ (86) 



— /xi / — 6/i3 1 — A; 2 1 + sn(x — / h(t)dt — at) 
M46= 2^ 4 cn(:r - J - at) 

I -6jU 3 1 - k~ cn(x - / h(t)dt - at) 

V /i 2 4 l + m(x-fh(t)dt-octy 1 j 



with 

Z -4^3 (l(l-k 2 )-l(l + k 2 ))-^ 
«=l ^ ^ ^ )• (88) 



— /Ji /— 6/i 3 k 2 1 + dn(x — / h(t)dt — at) 
UA1 ~ V /i 2 T ksn(x - J h(t)dt - at) 

I— 6/i3 fc 2 A;sn(:r — / h(t)dt — at) 
V /i 2 Tl + dn(a;-/^(t)d*-at)' ^ ^ 

with 

Q = / -^3(±3fc 2 + 2(l + ^ 2 ))-^ \ ^ (go) 



«48 = ~ ± 

2/i 2 



/-6/i 3 Ia2^' + dn 2 (x - / - at) 

V M2 ~ A; / -dn 2 (a;-//i(t)rft-at) 



, I S i , / ^ fc / -dn 2 (a;-/fe(t)dt-at) 
V ^ 1 J A/ + dn 2 (x -//»(*)* -a*)' 
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with 

, 4 W (±6/(i - mi + k' 2 ) + 2(1 + k'")) 
«=l 4^ )• ( 92 ) 



-/ii 24u 3 ,,A; / ± dn 2 (x - f /t(t)dt - at) 

M 49 = t; — =■= \ kl — = 

2^ 2 V ix 2 2 V / Fdn(x - / h(t)dt - at) 

, / -6/x 3 2^dn(a;- /fr(t)rft-at) 

± V + ^Fld^FTW^)' (93) 

with 



'-W 3 (±6 ) /±8*'(l + fc') 2 - (1 + A;' 2 ± 6tf)) - ^\ 
a=| — |. (94) 



— yUi —6/13, cn (x — f h(t)dt — at) + k' sn (x — f h(t)dt — at) 

M50= 2^ V + cn 2 (x - / fc(t)dt - at) - fc'sn 2 ^ - / h{t)dt - at) 



-6/x. 



cn 2 (x — / h(t)dt — at) — k'sn 2 (x — / h{t)dt — at) 



/i 2 cn 2 (:r — / h(t)dt — at) + /c'sn 2 (a; — / h(t)dt — at) 



with 

Q = / -4/i 2 /i 3 (±6(l-^)+2(l + fc0 2 )- / x 2 \ _ (g6) 



-IH -6fi 3 (l + kf 
W51 — — — ± 



2(l±fcsn(x-//i(t)dt-at)) 



2^2 V A*2 4 \\ (l + k)((l±sn(x-fh(t)dt-at))) 



(1 + jfc)((l ± sn(:r - / h(t)dt - at))) 



H 2 \\ 2(1 ±ksn(x-J h(t)dt- at)) 



(97) 



with 



-^(±12^(1 + fc) 2 + l + 6fc + P)- / u 2 
a=| — |. (98) 



^(1 + /csn(x - / /i(t)dt - at))(l + sn(x — / h{t)dt — at)) 
^ cn(x — / h(t)dt — at) 

± /~ 6//3 (l k)\j l ~ k cn{x - f h(t)dt - at) 

V A*2 V 2 ^(1 + ksn(x - f h(t)dt - at))(l + sn(:r - / h(t)dt - at)) ' 



with 

-/i 2 /i 3 (±6,/2(l - F)(l - k) 2 + 1 - Qk + A; 2 ) - /if 
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u 



2.5 



0.0 



10.0 



0.0 



Figure 1. Periodic wave solution for the mKdV equation which is given by family24 
(with (j,2 = 1) Ms = — 1) ~ 1 an d k = 0.5. 



Case 4: rational solution 




to x - f h(t)dt - (=g) t 



6/i 3 1 



(101) 



All the above solutions can be considered only in the case when /i 2 ^ 0. 

5. Applications to arterial mechanics 

In arterial mechanics j6], treating the arteries as a thin walled prestressed elastic tube 
with variable radius (or, with stenosis in [TJ) and blood as an inviscid fluid, the governing 
equation which models the weakly nonlinear waves in such a fluid-filled elastic tubes is 
the variable-coefficient mKdV equation flSJ). Also, in [i~8rf20] treating the arteries as 
a thin walled prestressed elastic tube and blood as an incompressible inviscid fluid , 
the governing equation which models the weakly nonlinear waves in such a fluid-filled 
elastic tubes is the mKdV equation ()5j). We obtained some useful solutions for these two 
equations such as periodic solution given by (157)) and solitary wave solution given by 
( l30l) and anti-kink wave solution given by ( l3Tl) which are shown in, (with h{t) = const.), 
fig{TJ figj2] and figJSJ respectively. 

6. Concluding remarks 

Using the classical symmetry analysis method we obtained two similarity transforma- 
tions (ITT)) and (l20l) that, for arbitrary h(t), transforms the original nonlinear PDE (JTJ 
into two nonlinear ODEs (ITS)) and (|2T)) respectively. Then we used the generalized Ja- 
cobi elliptic function expansion method to obtain many solutions for ()2TJ). The set of 
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Figure 2. Solitary wave solution for the mKdV equation which is given by family2 
(with fi2 — 1, (J-3 = 1 and h(t) = 1. 




Figure 3. Anti-kink wave solution for the mKdV equation which is given by family5 
(with fi2 — 1, (J-3 — — 1 and h{t) = 1. 

solutions from U\ to w 4 , from u 6 to itg, from « n to u 38 and -u 53 with h(t) = 0, can be 
found in [2TrBo] . To our best knowledge, the set of solutions u 5 , Ui and from u 39 to 
M52 are new solutions of (JT|) and are not shown in the current literature until now. 
The generalized variable-coefficient Gardner equation 

Ut + a{t)uu x + b(t)u 2 u x + hi{t)u xxx + d(t)u x + f(t)u = (102) 

was considered in (32] and solved with a new generalized algebraic method. Also, the 
set of solutions from u% to U4, from u G to u 9 , from iiu to ^33 and ^53 (with a(t) = fil, 
b(t) = fi2, h\(t) = fi3,d(t) = h(t) and f(t) = 0), are obtained in [32]. But, the set of 
solutions u 5 , u w and from u 39 to u 52 are not obtained in [32J. 
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Appendix A. Some special solutions of Equation ( I23j) . 



r p q (p 

1 e±C 

1 1 sinh(C) 
-1 1 cosh(C) 

1 -1 sin(C),cos(C) 

1 i 

1 -1 sech(C) 

1 1 csch(C) 

1 -2 1 tanh(C),coth(C) 
1 2 1 tan(C),cot(C) 

-11 sec(C),csc(C) 

1 -(1 + k 2 ) k 2 sn(C) 
k' 2 k 2 -k' 2 -k 2 cn(C) 

i i + k' 2 k' 2 ^(C) = ^g 

-k' 2 1 + k' 2 -1 dn(C) 

-k 2 k' 2 I -2k' 2 1 ^ 

k 2 -l 1+k 2 k 2 -l dn(C 



4 2 4 (l+fcsn(C)) 

1-k 2 1+k 2 1-k 2 cn(Q 

4 2 4 (l+sn(C)) 

1+k' 2 \£_ fcsn(C) 

4 2 4 (l+dn(C)) 

-*r -2{i + k' 2 ) (i+k-) |^ 

-4k' (l + k 2 ±6k') +(l + k') 2 2VEMO_ 



(fc'±dn^(0) 

:i + ^') 2 -2(1 + kr (i-kT ^jggjjg 

(l+fc) 2 l+6fc+fc 2 t, / (l+fc)(l=fcsn(C)) " 

4 4 ^ Y 2(l±fesn(C)) 

1-k 2 l-6fc+fc 2 (1-fc) 2 / l-fc cn(Q 

j v /(l+fcsn(C))(l+sn(Q) 



where fc 2 + k 
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